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HYPERSONIC  DUAL  COMBUSTION  RAMJET  ENGINES 


Richard  D.  Stockbridge,  Joseph  A.  Schetz,  Paul  J.  Waltrup, 
and  Frederick  S.  Billig 

The  Johns  Hopkins  University  Applied  Physics  Laboratory 
Laurel,  Maryland  20707 

INTRODUCTION 

Requirements  for  future  offensive  and  defensive  weapon  systems 

necessitate  the  development  of  long  range,  very  high  speed  missiles  to 
effectively  counter  the  continually  improving  capabilities  of  potentially 
hostile  nations.  Of  the  candidate  propulsion  cycles  available  to  power  these 
missiles,  only  rockets  and  advanced  ramjets  employing  supersonic  combustion  as 
their  primary  mode  of  combustion  are  capable  of  providing  the  hypersonic 
speeds  required.  However,  rockets  must  fly  exo-atmospheric  trajectories  to 

achieve  the  needed  ranges,  and  since  they  are  coasting  after  reentry,  their 
ability  to  make  corrections  and  to  intercept  maneuvering  targets  is  limited. 
On  the  other  hand,  advanced  hypersonic  ramjets,  which  remain  within  the 
atmosphere,  are  capable  of  sustained  powered  flight,  including  course  changes 
and  interception  of  maneuvering  targets. 

These  clear  advantages  of  an  advanced  hypersonic  airbreathing  pro¬ 
pulsion  system  establish  the  need  to  develop  a  basic  knowledge  and  under¬ 

standing  of  the  overall  engine  cycle,  individual  component  flowfields  and 
engine  thermochemistry  in  hypersonic  dual-combi  :ion  ramjet  engines.  This 

goal  is  being  pursued  by  both  an  analytical  approach  and  by  an  experimental 

approach  that  focus  on  different  aspects  of  this  problem.  The  analytical  ap¬ 
proach  involves  the  development  of  theoretical  models  that  describe  the  dual¬ 
combustion  process  with  empnasis  on  the  supersonic  combustor  and  exit 

nozzle.  The  experimental  approach  is  being  used  to  establish  a  sound  data 

base  and  a  corresponding  semi-empirical  model  applicable  over  the  range  of 
simulated  flight  conditions  that  will  describe  the  interaction  of  the 

compression  field  generated  at  the  entrance  of  supersonic  combustors  with  the 
exit  flowfield  of  hypersonic  air  inlets. 


Dump 

M  1  M  >  1  combustor 


Fig.  1  Schematic  of  dual  combustion  ramjet  engine. 


Fig.  1  is  a  schematic  diagram  of  the  Dual  Combustion  Ramjet  (DCR) 
engine  concept.  The  freestream  air  is  initially  compressed  by  an  external  in¬ 
let  diffuser.  The  inner  cowl  lip  then  subdivides  the  compressed  air,  and  a 
small  portion  is  ducted  to  a  small,  dump-type  subsonic  combustor.  All  of  the 
fuel  is  added  in  the  subsonic  combustor,  which  acts  as  a  fuel-rich,  hot-gas 
generator  for  the  supersonic  combustor.  The  major  portion  of  the  air  bypasses 
the  gas  generator  and  is  ducted  supersonically  through  an  air  conduit  to  the 
supersonic  combustor,  where  it  mixes  and  burns  with  the  exhaust  of  the  gas 
generator.  The  final  shocks  shown  in  the  supersonic  duct  (depicted  as  two 
sets  of  intersecting  oblique  waves)  result  from  the  combustion-induced  pres¬ 
sure  disturbances  generated  by  the  supersonic  combustion  process.  The  region 
is  commonly  referred  to  as  the  combustor/ inlet  interaction  region. 


COMBUSTOR/ INLET  INTERACTIONS 


Experimental  Program 


Data  from  connected-pipe  tests  of  Dual  Combustion  Ramjet  combustors 
h  ve  shown  the  existence  of  a  precorabus tion  pressure  rise  beginning  in  the 
annular  supersonic  flow  upstream  of  the  gas  generator  base  and  extending  into 
the  supersonic  combustor  (Ref.  1).  This  precombustion  pressure  rise  is  being 
simulated  in  a  simplified  test  setup  using  overexpanded  cold  air  in  order  to 


develop  a  data  base  over  a  wide  range  of  test  conditions  and  geometries. 
Fig.  2  on  p.  7  Is  a  schematic  and  photograph  of  the  DCR  Combus tor/ Inlet 
Interaction  (CII)  test  hardware.  The  annular  isolator  duct  that  conducts 
supersonic  air  to  the  supersonic  combustor  surrounds  the  sonic  or  supersonic 
exhaust  flow  from  the  gas  generator  simulator.  Test  variables  include  the 
annulus  Mach  number  (1.75,  2.4  and  3.0)  and  Reynolds  number,  the  gas  generator 
exit  Mach  number  (1.0  to  2.5),  the  airflow  split  between  the  supersonic 
annulus  and  gas  generator  (1:1  to  6:1),  and  the  length  of  constant  area  duct 
downstream  of  the  exit  of  the  gas  generator.  The  primary  instrumentation 
includes  wall  and  instream  pressure  measurements. 

The  data  from  this  experimental  program  will  be  used  to: 

(1)  yield  a  data  correlating  equation  that  will  define  the  location 
and  extent  of  the  precombustion  compression  field  in  the  annu¬ 
lar  isolator  duct  that  conducts  compressed  air  from  the  engine 
inlet  to  the  supersonic  combustor  as  a  function  of  the 
appropriate  variables; 

(2)  describe  the  boundary  layer  growth,  statistics  of  the  turbu¬ 
lence  and  average  fiowfield  property  radial  profiles  in  the 
annular  isolator  duct;  and 

(3)  measure  pitot  pressure  radial  profiles  in  the  supersonic 
combustor  section  at  various  axial  stations  downstream  of  the 
gas  generator  base. 

Data  correlations  from  the  CII  test  program  will  be  used  in  conjunction  with 
the  DCR  combustor  test  data  to  design  inlet  air  ducts  which  will  prevent  the 
precorabustion  pressure  disturbance  from  disrupting  the  external  air  inlet 
fiowfield. 

Resu 1 ts 

Various  procedures  have  been  evaluated  to  calculate  the  instream 
flow  properties  in  a  supersonic  annular  fiowfield  from  measured  instreara  pitot 


pressures,  wall  pressures  and  plenum  conditions.  Each  procedure  achieves 
analytical  closure  on  all  flowfield  property  variations  by  deriving  the  values 
for  a  second  instream  flowfield  property  as  a  function  of  the  radial  (r)  and 
axial  (2)  position  based  on  various  theoretical  premises.  One  procedure 
assumes  a  constant  static  pressure  across  the  annular  duct  equal  to  the  aver¬ 
age  of  the  five  measured  wall  pressures  at  that  axial  station.  Another  pro¬ 
cedure  assumes  a  constant  static  pressure  which  conserves  mass  at  that 
station.  A  third  variation  assumes  a  static  pressure  profile  at  each  axial 
position  that  is  quadratic  with  respect  to  the  radius  r,  conserves  mass  and 
matches  the  measured  static  pressures  at  the  wall.  All  of  these  procedures 
were  found  to  be  unsatisfactory  (Refs.  2,  3).  One  reason  for  this  was  that 
some  of  the  computed  instream  total  pressures  (pt)  at  various  radii  were 
higher  than  the  measured  plenum  total  pressure,  thereby  violating  entropy 
constraints.  Two  additional  unsatisfactory  results  of  these  initial  assump¬ 
tions  were:  (1)  the  calculated  displacement  thicknesses  and  turbulence 

intensities  were  unrealistic;  and  (2)  the  correlation  of  the  fluctuating 
components  of  density  and  velocity  (avg.  of  p  u)  close  to  the  walls  was 
greater  than  zero.  These  difficulties  were  exacerbated  by  the  high 

sensitivity  of  Mach  number  and  total  pressure  to  small  variations  in  the 
assumed  instream  static  pressure  profile. 

A  new  procedure  has  been  developed  that  (1)  conserves  mass;  (2) 
matches  the  measured  wall  pressure;  (3)  ensures  that  the  stagnation  entropy 
never  decreases  in  the  flow  direction;  and  (4)  correctly  models  the  effects  of 
viscosity  both  in  the  center  of  the  duct  and  next  to  the  walls.  This  proce¬ 
dure  assumes  a  symmetric  radial  total  pressure  profile,  pt(r,  z),  that  is 

constant  and  equal  to  the  plenum  total  pressure  (p.  )  over  a  central  inviscid 

co 

core  at  the  center  of  the  duct,  and  is  quadratic  with  respect  to  distance  from 
the  midpoint  between  the  annular  walls  for  the  two  wall  regions  (see  Fig.  3  on 
p.  9).  The  quadratic  curve  fit  allows  three  constraints  to  be  imposed:  (1) 
zero  slope  at  the  junction  of  the  quadratic  with  the  central  inviscid  core 
segment,  (2)  matching  of  the  centerline  total  pressure  (pt  ),  and  (3)  matching 
of  the  wall  pressure.  The  position  of  the  junction  point  is  determined 
iteratively  to  conserve  mass.  This  new  integral  analysis  overcomes  all  of  the 
deficiencies  of  the  previous  methods  while  retaining  their  desirable  proper¬ 
ties  (conserving  mass  and  matching  the  measured  wall  static  pressures). 
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Radial  position  from  center  of  annulus  to  either  wall 


Assumed  total  pressure  radial  profiles  in  the  Cl  I  annular  duct  normalized  with 
respect  to  the  plenum  total  pressure  calculated  to  match  the  plenum  and  wall 
pressures  and  to  conserve  mass. 


Fig.  3  on  page  9  shows  pt(r,  z)  for  the  three  axial  positions  (z)  in  the 
annular  isolator  duct  for  which  pitot  pressure  p^  measurements  were  made. 
The  variable  z  is  the  distance  from  the  annular  nozzle  throat  to  the 
designated  axial  position.  Fig.  3  on  page  9  shows  graphically  the  growth  of 
the  viscous  wall  region  and  the  decay  of  che  inviscid  core  as  the  flow 
progresses  downstream. 

The  determination  of  the  boundary  layer  thickness,  6,  displacement 
thickness,  6*,  momentum  thickness,  0,  and  kinetic  energy  thickness,  6  , 

affords  a  quantitative  measurement  of  the  importance  of  viscosity  in  the  DCR 
combustor/ inlet  isolator  duct.  These  data  are  important  in  the  determination 
of  the  required  geometric  area  divergence  for  the  Isolator  duct  in  the  DCR 
engine  and  in  the  development  of  a  da ta-correla ting  equation  that  specifies 
the  location  of  the  precombustion  compression  field  in  the  isolator  duct  as  a 
function  of  various  parameters  that  may  include  6*  and/or  0  (Ref.  4).  The 
equations  defined  in  Ref.  5  for  these  thicknesses  are  listed  in  Table  1  on  p. 
11.  They  are  for  compressible,  axisymmetric  flow  around  or  within  a  body  with 
transverse  radius  R.  For  the  annulu3  inner  wall  boundary  layer  thicknesses, 
the  flowfield  is  external  and  R  is  positive,  whereas  the  opposite  applies  for 
the  annulus  outer  wall. 

Fig.  4  on  p.  12  shows  the  calculated  boundary  layer  thickness 
plotted  against  z.  There  was  no  change  in  these  boundary  layer  results  when 
the  Reynolds  number,  Re,  was  varied  from  12.1  to  19.0  x  107/m.  At  the  end  of 
the  0.513  m  long  constant  area  annular  duct  (z  ■  0.607  m),  the  fractions  of 

A 

the  total  cross-sectional  area  (0.016  nr)  occupied  by  the  various  boundary 
layer  thicknesses  are  74%  for  6,  13%  for  5*,  7.6%  for  and  4.2%  for  0. 

As  shown  in  Figs.  5-8  on  p.  13  and  14,  this  boundary  layer  growth 
affects  the  mean  flowfield  property  variations  with  respect  to  z.  At  the 
center  of  the  duct,  for  z  increasing  from  0.275  to  0.506  m  downstream  of  the 
nozzle  throat,  the  Mach  number,  M,  falls  from  2.42  to  2.33,  the  velocity,  u, 
falls  from  565  to  555  m/sec  and  the  mass  flux,  pu  (p  ■  density),  rises  from 
1220  to  1308  kg/ra4-sec.  Mote  also  the  effect  of  the  boundary  layer  growth  on 
the  average  values  of  these  properties  in  the  boundary  layer.  As  expected, 
both  the  Mach  number  and  velocity  are  retarded  at  a  fixed  distance  from  either 


Table  1 

Defining  equations  for  integral  thicknesses. 


Displacement  thickness  5*: 
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1  +  R  dV 


Momentum  defect  thickness  d: 


Kinetic  energy  defect  thickness  S« E 


transverse  body  radius  whose  sign  is  positive  for  external 

and  negative  for  internal  flows 
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Fig.  6  Average  velocity  vs.  radial  position  in  the  Cl  I  Mach  2.4  annular  flowfield 


Normalized  radius 


Average  Pitot  pressure  vs.  radial  position  in  the  Cl  I  Mach  2.4  annular  flowfield. 


wall  as  the  boundary  layer  grows  in  the  downstream  direction.  The  non-smooth 
pattern  for  pu  in  the  boundary  layer  is  also  evident  in  the  average  pitot 
pressure  data  (p^  ,  Fig.  8  on  p.  14),  which  is  the  experimentally  measured 
data  upon  which  all  property  calculations  depend.  While  the  p^  curve  for  z 
*  0.275  m  shows  a  continuous  rise  in  p^  with  distance  from  either  wall,  the 
departure  from  this  expected  behavior  is  moderate  for  z  **  0.398  m  and  more 
severe  at  z  =*  0.506  m  in  bo th  wall  regions.  Thus,  the  consistency  of  this 
pattern  with  respect  to  variations  in  z  and  for  either  the  inner  or  outer  wall 
supports  the  validity  of  the  measured  results. 

The  results  for  M  and  u  in  Figs.  5  and  6  on  p.  13  are  invariant  with 
respect  to  the  above  Reynolds  number  range.  The  pu  and  p^  plots  in  Fig.  7  and 
8  on  p.  14  are  for  Re/m  =*  12.1  x  107  and  a  plenum  total  pressure  pt  *  1241 
kPa.  The  curves  for  these  quantities  at  Re/m  *  19  x  107  are  similar  in  every 
respect  except  absolute  magnitude. 

In  order  to  quantify  the  magnitude  of  the  fluctuating  component  of 
various  flowfield  properties  in  the  CII  annular  isolator  duct,  the  root  mean 
squares  (rms)  of  the  time  series  for  the  fluctuating  component  of  the  pitot 
pressure  p^,  velocity  u,  density  p  and  mass  flux  pu  have  been  computed.  Note 
that  the  time  average  of  the  fluctuating  component  of  any  property  time  series 
is  zero.  These  computed  results  have  been  normalized  with  respect  to  the 

average  value  of  the  specified  property  at  the  radial  position  in  the  annular 

duct  for  which  the  velocity  is  a  maximum  for  that  axial  station.  Figs.  9-12 
on  p.  16-17  show  these  data,  respectively,  plotted  against  the  normalized 
radial  position  where  the  data  were  measured  for  three  axial  positions  in  the 
CII  isolator  duct,  z  m  0.275,  0.398  and  0.506  m  downstream  of  the  annular 
nozzle  throat.  The  last  station,  z  *  0.506,  is  0.1  m  upstream  of  the  gas 

generator  base.  In  these  figures,  a  fluctuating  component  of  a  flowfield 
property  is  designated  by  a  superscript  tilda  and  an  average  value  by  a 
superscript  bar. 

These  results  show  that  the  intensity  of  the  turbulence  is  highest 
in  the  boundary  layer  and  increases  as  the  boundary  layer  thickens.  In  the 

CII  isolator  duct,  viscosity  and  its  various  ramifications,  such  as  signif¬ 
icant  boundary  layer  growth,  boundary  layer  influence  on  the  mean  flow  and 
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Fig.  9  Pitot  pressure  fluctuations  in  the  Mach  2.4  ducted  flowfield. 
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Fig.  10  Velocity  fluctuations  in  the  Mach  2.4  ducted  flowfield. 
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Fig.  1 1  Density  fluctuations  in  the  Mach  2.4  ducted  f lowfield. 


Fig.  12  Mass  flux  fluctuations  in  the  Mach  2.4  ducted  f lowfield. 


Fig.  13  shows  that  the  statistical  correlation  of  p  and  u  in  the 
duct  is  zero  outside  of  the  boundary  layer  and  negative  inside  the  boundary 
layer.  This  is  an  expected  result,  typical  of  turbulent  boundary  layers,  and 
while  primarily  of  theoretical  interest,  does  offer  an  important  confirmation 
of  the  validity  of  the  computational  procedure. 

The  investigation  of  the  simulated  precombustion  compression  field 
in  the  overexpanded  annular  flowfield  is  the  primary  purpose  of  this  experi¬ 
mental  investigation.  The  calculation  of  5,  <S*,  0  and  6,,_  throughout  the 
annular  duct  from  measured  data  has  been  completed  primarily  so  that  one  or 
more  of  these  thicknesses  can  be  incorporated  into  the  correlating  equation 
that  will  be  derived  to  predict  the  precorabus tion  compression  field 
location.  Ref.  4  shows  that  0  was  an  Important  parameter  in  an  analogous 


correlation,  for  the  pare  scrarajet  geometry.  The  location  and  length  of  the 

simulated  precombustion  compression  field  is  a  function  of  t  e  following  test 

variables:  (1)  the  degree  of  overexpansion  relative  to  ambient  exhaust 

conditions;  (2)  the  mass  flow  ratio  (w  /w  )  between  the  outer  annular 

a  g 

and  inner  gas  generator  flows;  (3)  the  gas  generator  exhaust  nozzle  Mach 

number  and  contour;  (4)  the  length  of  duct  that  extends  downstream  from  the 

end  of  the  annulus;  and  (5)  the  annular  nozzle  Mach  number.  These  items 

are  being  systematically  varied  over  the  entire  experimental  range  of 

interest.  Figs.  14-16  on  p.  20-22  show  various  measured  wall  pressure 

distributions  that  simulate  the  precombustion  compression  field.  A  pressure 

rise  ratio  greater  than  the  normal  shock  pressure  rise  ratio  for  the  actual 

nozzle  exit  Mach  number  (2.4)  has  been  obtained  (Fig.  14  on  p.  20)  for  a  case 

where  w^/w^  =  0.506.  This  actual  mass  flow  split  corresponds  to  an  effective 

mass  flow  split  of  1.04  when  the  effect  of  the  maximum  temperature  difference 

between  the  two  streams  in  the  operating  engine  is  included.  Fig.  16  on  p.  22 

shows  the  isolated  effect  of  varying  w  /w  .  For  high  w  ,  the  gas  generator 

a  g  g 

flow  expands  supersonically,  lowers  the  base  pressure  and  has  the  aspirating 
effect  of  translating  downstream  the  position  of  the  compression  field  in  the 
annulus.  This  allows  a  greater  degree  of  overexpansion  and,  consequently,  a 
higher  overall  pressure  rise  ratio. 


Discussion 

The  value  of  the  boundary  layer  integral  thicknesses  6*,  0, 

and  ($Kg  are  very  sensitive  to  the  value  of  the  boundary  layer  thickness  6, 
which  is  usually  defined  to  be  the  distance  from  the  wall  at  which  the  local 
average  velocity  is  equal  to  0.99  times  its  "f rees tream"  value.  This  sensi¬ 
tivity  enters  the  equations  in  Table  1  on  p.  11  through  the 

and  u^  quantities  in  the  integral  equations  and  through  the  6  limit  of 
integration.  For  an  internal  flowfield  with  a  static  pressure  gradient  normal 
to  the  wall  (the  CII  flowfield),  there  is  not,  in  general,  a  definitive  value 
for  6  because  there  is  no  definitive  "freestream"  velocity  (Ref.  5,  p.  19). 
Because  the  integral  thicknesses  are  very  sensitive  to  the  value  of  6,  it  is 
recommended  in  Ref.  5  for  supersonic  flows  that  6  be  defined  as  the  thickness 
of  the  region  next  to  the  wall  in  which  the  local  total  pressure  increases  to 
99%  of  its  original  (plenum)  value.  This  definition  removes  the  ambiguity 
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section  downstream  of  the  annulus. 


that  would  arise  if  the  velocity  ratio  were  used  and  provides  a  common  basis 
of  comparison.  It  is  also  philosophically  consistent  with  the  total  pressure 
quadratic  curve  fit  analytical  closure  procedure  that  is  used  here  to 
calculate  instream  properties.  On  the  basis  of  these  arguments  and  on  the 
numerical  results  of  calculating  the  iategral  thicknesses  based  on  different 
criteria  for  the  selection  of  6,  this  total  pressure-based  5  definition  has 
been  used  consistently  in  the  integral  thickness  calculations  presented  in 
this  report. 

The  CII  experimental  results  are  compared  in  Tables  2-4  on  p.  24-25 
to  compressible  turbulent  boundary  layer  data  obtained  by  other  inves¬ 
tigators.  In  these  tables,  y  is  the  distance  from  the  wall.  When  used  as  a 
subscript,  5  indicates  the  value  of  a  property  at  the  boundary  layer  edge. 
All  of  the  turbulent  fluctuation  data  in  the  surveyed  literature  (Refs.  6-9) 
came  from  velocity  measurements  using  hot  wire  anemometers.  The  comparable 
CII  data  came  from  measuring  the  pitot  pressure  with  both  a  high  response 
Kulite  brand  miniature  pressure  transducer  and  with  Alinco  brand  pressure 
transducers.  Although  a  direct  comparison  of  the  flow  property  statistics 
determined  from  a  hot  wire  and  from  these  pressure  transducers  for  the  same 
flowfield  has  not  been  made,  the  CII  data  in  Tables  2-4  on  p.  24-25  compares 
favorably  with  the  other  data,  particularly  in  the  supersonic  Mach  number 
range  for  the  velocity  data.  Table  5  on  p.  26  summarizes  some  of  the  boundary 
layer  data  at  the  3  axial  locations  in  the  CII  annular  duct  where  pitot 
pressure  radial  profiles  were  measured.  The  Reynolds  number  Re^  range  was  40- 
68  x  10"3. 

The  Kulite  brand  miniature  pressure  transducer  frequency  response  is 
flat  to  20  kHz,  whereas  the  cutoff  frequency  of  the  Alinco  is  only  75  Hz. 
Both  are  silicone-diaphragm,  strain-gauge  transducers.  The  Kulite  transducer 
has  no  internal  volume  (the  diaphragm  is  an  external  surface  of  the 
transducer);  it  was  mounted  in  the  flow  at  the  tip  of  a  pitot  probe.  The 
Alinco  transducers  have  0.25  in  Internal  volume  between  the  entrance  port  and 
the  pressure-sensitive  diaphragm;  they  were  connected  to  the  pitot  probes  with 
about  1  ft  of  0.065  in.  inner  diameter  Nylon  tubing,  which  added  an  additional 
0.04  in  or  16%  to  the  internal  volume.  The  Kulite  was  positioned  in  the 
approximate  center  of  the  inner  boundary  layer  at  a  normalized  radial 


Comparison  of  turbulent  density  fluctuation  p'  results  with  other  compressible 
turbulent  boundary  layer  data. 


m5 

rms  p'  (y)/ 

Avg  p  (y) 

Max*  rms  p'  (y )/ 

Avg  p  (5) 

(Z,  cm) 

Max* 

(y  IS) 

at  y  =  5 
(5,  mm) 

(y/6) 

Stockbridge 

2.38 

0.090 

0.088 

(27.5) 

(0.217) 

(0.217) 

2.37 

0.147 

0.140 

(39.8) 

(0.156) 

(7.3) 

(0.156) 

2.33 

0.173 

0.023 

0.146 

Owen 

(50.6) 

(0.370) 

(8.9) 

(0.217) 

7.2 

0.10 

(0.43) 

0.067 

(33) 

0.062 

(0.56) 

Laderman 

9.4 

0.137 

(0.543) 

0.026 

(141) 

*Max  refers  only  to  the  maximum  of  the  available  tabular  data. 
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Table  4 

Comparison  of  (w/A)'  results  with  other  compressible  turbulent 
boundary  layer  data. 


m5 

(Z,  cm) 

rms  (w/A)'  (y)/ 

Avg  (w/A)  (y) 

Max*  rms  (w/A)'  (y)/' 

Avg  (w/A)  (5) 

Max* 

(y/5) 

at  y  =  5 
(5,  mm) 

(y/5) 

Stockbridge 

2.38 

0.068 

0.012 

0.060 

(27.5) 

(0.541) 

(5.3) 

(0.541) 

2.37 

0.086 

0.014 

0.069 

(39.8) 

(0.389) 

(7.3) 

(0.389) 

2.33 

0.116 

0.019 

0.088 

(50.6) 

(0.320) 

(8.9) 

(0.320) 

Owen 

7.2 

0.15 

0.065 

0.025 

(0.04) 

(33) 

(0.04) 

*Max  refers  only  to  the  maximum  of  the  available  tabular  data. 
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Table  5 

Cl  I  boundary  layer  data  at  three  axial  positions  in  the  annular  isolator  duct. 


Distance  from  nozzle  throat,  z  [m] 


0.275 

0.398 

0. 

5  [mm] 

5.3 

7.3 

8. 

5*  [mm] 

1.0 

1.1 

1. 

d  [mm] 

0.3 

0.5 

0. 

5ke  [mm] 

0.6 

0.9 

1. 

m6 

2.38 

2.37 

2. 

Hee  X  1CT3 

40.4 

61.8 

68. 

T5  [°K] 

139 

139 

139 

rms  u'(5)  [m/sec] 

1.58 

1.92 

2. 

rms  p'(5)  [kg/m3] 

0.032 

0.032 

OJ 

rms  (pu)'  (5)  [kg/m2  sec] 

15.29 

17.62 

23. 

u  (5)  [m/sec] 

560.5 

559.6 

554. 

p  (S)  [kg/m3] 

2.194 

2.178 

2. 

Note:  All  values  for  integral  boundary  layer  thicknesses  are  the  largest  for 
both  walls. 


position,  R,  of  -0.79  and  accounts  for  the  spike  in  the  property  fluctuation 
plots  at  that  radius.  It  is  significant  that  the  maximum  values  for  rms  u' 
and  rms  p'  were  measured  with  Alinco  transducers  at  radii  closer  to  both 
walls  than  R  *  -0.79  because  it  shows  that  the  limited  bandwidth  of  the  Alinco 
transducers  did  not  change  the  trend  of  the  correlation  of  turbulence 
intensity  with  proximity  to  either  annular  wall. 

The  Kistler  data  (Ref.  6)  in  Table  2  on  p.  24  were  obtained  at  a 
single  station  on  the  flat  wall  opposite  flexible  two-dimensional  (2D)  nozzle 
plates.  The  Reynolds  number  ReQ  range  was  26  to  35  x  10  .  The  static 
pressure  was  constant  through  the  boundary  layer,  and  the  test  section  walls 
were  adiabatic. 

The  Waltrup  and  Schetz  data  (Ref.  7)  in  Table  2  on  p.  24  is  for  a  2D 
turbulent  boundary  layer  along  a  test  section  wall  that  is  subjected  to  an 
adverse  pressure  gradient  by  the  insertion  of  a  ramp  into  a  2D  nozzle.  The 

Reynolds  number  Re^  range  was  20  to  40  x  10  and  the  walls  were  adiabatic. 

The  Owen  data  (Ref.  8)  in  Tables  2,  3,  and  4  on  p.  24-25  came  from 
the  Mach  7.2  external  flowfield  about  an  axisymmetric  body.  The  test  boundary 
layer  was  formed  on  a  10°  semiapex  angle  cone-ogive-cylinder  3.3  m  long  and 
0.203  ra  in  diameter.  The  wall  temperature  was  only  half  the  recovery 

_3 

temperature  and  Re^  ■  13  x  10  . 

The  Laderraan  data  (Ref.  9)  in  Tables  2  and  3  on  p.  24-25  came  from  a 
boundary  layer  that  was  formed  on  the  upper  wall  of  a  contoured  2D  flexible 
nozzle  and  test  section.  There  was  severe  heat  transfer  to  the  wall;  the  wall 
temperature  was  0.4  times  the  recovery  temperature.  The  test  surface  was 

polished  to  a  micro-in.  finish  and  the  wall  was  actively  cooled.  The 

Reynolds  number  Re^  was  40  x  10  ,  and  the  Mach  number  was  9.4. 

COMPONENT  AND  ENGINE  CYCLE  MODELING 

Rather  than  develop  a  single  model  describing  the  complete  DCR 
supersonic  combustor  and  exit  nozzle  flow  processes,  one  which  would  be 
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extremely  complex  and  prohibitively  expensive  to  use,  an  alternative, 
"building  block"  approach  has  been  adopted.  This  approach  permits  each  dis¬ 
tinct  region  (Fig.  17)  downstream  of  the  combustor  entrance,  viz.,  the  base 
flow,  central  core  flow  and  wall  boundary  layer,  to  be  described  separately 
and  then  matched  to  each  other.  Since  each  region  is  amenable  to  certain 
simplifying  physical  assumptions  not  necessarily  applicable  in  another  region, 
this  approach  permits  a  simpler,  yet  physically  realistic  description  of  each 
region  compared  to  an  all-encompassing  analysis,  and  a  complete  solution  is 
orders  of  magnitude  less  expensive  to  obtain. 

Fig.  17  illustrates  the  region  of  interest  in  these  analyses.  The 
shocks  shown  in  the  supersonic  duct  result  from  the  combustion-induced  pres¬ 
sure  disturbances  generated  by  the  mixing  and  heat  release  processes.  The 
region  is  commonly  referred  to  as  the  corabus tor/inlet  interaction  region. 
With  no  heat  release,  the  shock  structure  is  not  present  and  expansion  waves 
are  possible  depending  on  the  geometry  and  flow  conditions  in  the  inner  jet. 
With  heat  release,  the  pressure  rise  associated  with  this  shock  system  can 
vary  from  low  values  equivalent  to  those  corresponding  to  oblique  waves  up  to 
the  value  corresponding  to  a  normal  shock  in  the  incoming  air.  Thus,  at  low 
'iach  number  and  high  heat  release  rates,  the  flow  in  the  laain  combustor  is 
subsonic,  which  is  the  mode  of  operation  in  the  conventional  "scrarajet"  or 
"dual-mode"  engine  (Ref.  10). 


A  modular  approach  to  the  analysis  of  this  complex  fiowfield  has 
been  adopted.  Fig.  17  on  p.  28  also  shows  the  regions  of  flow  that  are 
handled  in  these  modules.  The  mixing  and  burning  of  the  central  portion  of 
the  flow  is  the  main  module  and  will  ultimately  use  input  from  modules  that 
treat  the  shock  structure  in  the  incoming  outer  flow  and  from  a  module  that 
handles  a  detailed  treatment  of  the  "base  region"  of  the  injector.  A  separate 
module  is  used  to  calculate  the  wall  boundary  layer  and  then  is  coupled  to  the 
main  flow  solution.  This  approach  is  made  necessary  by  the  substantial 
difference  in  grid  resolution  required  in  these  two  parts  of  the  flow.  A  much 
finer  grid  is  required  in  the  boundary  layer  than  in  the  mixing  and  burning 
zone,  so  that  wall  properties  can  be  accurately  predicted. 


One  of  the  important  complicating  features  of  the  flow  in  the  com¬ 
bustor  is  the  recirculation  region  in  the  base  of  the  thick  "lip"  of  the  gas 
generator  exhaust  (Fig.  17  on  p.  28).  The  flow  in  this  region  interacts  with 
the  hot  exhaust  jet  and  the  inlet  air  to  determine  the  pressure  and  the  other 
flow  variables  at  the  beginning  of  the  mixing  and  burning  region  in  the  main 
combustor.  The  accurate  prediction  of  the  subsequent  flow  in  the  combustor  is 
clearly  dependent  upon  the  determination  of  good  "initial  conditions"  as  a 
result  of  an  analysis  of  this  base  flow. 


Fig.  18  Schematic  illustration  of  base  flow  entrainment  model. 
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The  flowfield  for  a  simple,  planner  base  is  shown  schematically  in 
Fig.  18  on  p.  29.  The  analysis  of  this  complex  flow  has  been  compart¬ 

mentalized.  The  flow  in  the  upstream  boundary  layer  and  some  of  the  inviscid 
flow  above  the  boundary  layer  is  taken  inviscidly  through  the  corner  expansion 
and  the  lip  shock.  None  but  the  most  complicated  existing  analyses  include 

the  lip  shock.  The  subsequent  downstream  flow  outside  of  the  viscous  base 
region  is  also  treated  as  inviscid.  The  viscous  treatment  of  the  flow  in  the 
base  begins  at  the  viscous  throat  and  proceeds  upstream.  The  flow  from  the 
viscous  throat  to  the  rear  stagnation  point  is  analyzed  using  generalized 
versions  of  the  methods  in  Refs.  11-14.  The  flow  in  the  region  above  the 

dividing  streamline  through  the  rear  stagnation  point  and  upstream  from  the 
rear  stagnation  point  to  the  downstream  end  of  the  constant  pressure  region  is 
treated  using  integrated  equations  of  motion  but  with  the  pressure  gradient 
taken  as  predicted  for  a  centered  wave  pattern.  A  unique  composite  solution 
for  the  whole  is  selected  by  requiring  continuity  of  the  mass  flow  in  the 
shear  layer  from  the  downstream  and  upstream  proceeding  portions  of  the 

ana lysis . 


Some  comparisons  of  predictions  and  experiments  for  the  planar  base 
flow  are  given  in  Figs.  19  to  21  on  p.  31-33.  Figs.  19  and  21  show  predicted 
and  experimental  centerline  pressure  distributions  at  two  upstream  Mach 
numbers  in  the  range  of  interest.  Clearly,  the  base  pressure  and  the  major 
features  of  the  flow  are  well  predicted.  Fig.  20  on  p.  32  shows  that  the 
important  effects  of  the  size  of  the  upstream  boundary  layer  are  adequately 
accounted  for. 

The  main  mixing  and  burning  calculations  can  currently  be  made  with 
either  a  greatly  extended  version  of  the  code  presented  in  Ref.  15  or  the  Los 
Alamos  VNAP-II  code  developed  by  M.  C.  Cline.  The  former  is  based  on  the 
boundary  layer  (parabolic)  formulation  of  the  equations  of  motion  and  the 
latter  is  a  full  axlsymmetrtc  solution  of  the  Navter-S tokes  equations,  which 
permits  the  Important  effects  of  pressure  feedback  and  radial  pressure  gradi¬ 
ents  to  be  included.  The  extensions,  in  each  case,  include  the  integrated, 
turbulent  kinetic  energy  (TKE)  method  for  calculating  the  eddy  viscosity 
(following  Ref.  lb),  the  doub le - f lameshee t  model  of  Ref.  17,  the  chemical 
equilibrium  calculation  procedure  of  Ref.  13,  and  the  integral,  confined  flow 
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Fig.  20  Effect  of  upstream  boundary  layer  momentum  thickness  on  base  pressure. 


Fig.  21  Predictions  and  experiments  for  pressure  distribution  behind  a  base  at  M1  =  2.30. 


with  heat  release  treatment  of  Refs.  10  and  19.  For  the  parabolic  formula¬ 
tion,  the  first  extension  was  to  accommodate  confined  mixing  in  a  duct  with  a 
central  fuel  jet  rather  than  jet  mixing  in  an  infinite  stream  as  had  been  done 
in  Ref.  20.  This  was  accomplished  by  holding  the  stream  function  ij;  along  the 

outer  boundary  fixed  at  its  initial  value  in  order  to  maintain  a  constant 

iIAa 

total  mass  flow  in  the  duct.  Velocity  "slip"  was  allowed  along  the  duct 
wall.  Two  options  are  incorporated:  in  option  1,  p(z)  is  specified  and  the 
cross-sectional  area  A(z)  is  computed;  in  option  2,  the  inverse  is  done.  The 
second  extension  incorporated  a  more  advanced  turbulent  transport  model  in 
both  of  these  options.  An  integrated  equation  for  turbulent  kinetic  energy 
was  used  to  predict  the  variation  of  eddy  viscosity  along  the  duct  following 
the  method  in  Ref.  16.  The  third  and  largest  extension  incorporated  equilib¬ 
rium  chemical  reactions  into  the  mixing  zone. 

The  C-H-O-N  chemical  system  is  written  to  include  the  following 
species:  C>2,  N2,  H2 ,  CO,  CO2,  OH,  0,  H,  N,  and  NO.  Because  of  the  high 

static  pressures  and  temperatures  and  the  need  to  reduce  the  problem  to 
manageable  proportions,  local  diffusion-controlled  chemical  equilibrium  was 
assumed.  The  computer  code  given  in  Ref.  18  was  modified  extensively  to 
reduce  the  cost  per  calculation  for  an  enthalpy-pressure  entry  by  an  order  of 
magnitude,  to  about  one  cent  per  point.  Still  further  simplifications  were 
necessary  to  bring  the  computation  cost  to  an  acceptable  level.  For  example, 
about  25  radial  points  with  1000  axial  stations  are  required  to  obtain 
acceptably  accurate  solutions  in  a  0.1-m  radius,  2-m-long  duct.  At  one  cent 
per  point,  the  cost  of  the  equilibrium  calculations  is  about  $250-$500,  com¬ 
pared  to  less  than  $50  for  the  remainder  of  the  entire  mixing  calculations. 

To  decrease  the  number  of  required  equilibrium  calculations,  the 
double-f lamesheet  model  of  Ref.  17  was  adapted  to  the  present  problem.  Simply 
stated,  the  model  comprises  an  inner  H2  flame  and  an  outer  CO  flame  with  only 
the  O2,  H2,  N2,  H2O,  CO,  and  CO2  3pecies  assumed  to  be  present.  As  long  as 
the  temperature  remains  below  about  2500oK,  the  model  provides  a  good  approxi¬ 
mation  to  the  exact  equilibrium  state.  Therefore,  a  composite  model  was 
adopted.  At  each  axial  station,  a  complete  calculation  is  made  using  the 
double-f lamesheet  model.  The  temperature  profile  is  searched  for  points 
having  a  predicted  temperature  greater  than  2500°K.  For  those  points,  a 


complete  equilibrium  calculation  is  made  before  moving  to  the  next  axial 
station.  The  resulting  cost  reduction  factor  has  been  about  three  to  five. 

The  fourth  extension  was  necessary  because  the  hot  combustion  prod¬ 
ucts  from  the  main  mixing  and  burning  region  are  in  the  shear  layer  at  the 
outer  edge  of  the  wall  boundary  layer,  and  they  significantly  affect  wall  heat 
transfer  and  3kin  friction.  Furthermore,  the  turbulence  produced  in  that 
shear  layer  influences  the  wall  boundary  layer,  just  as  high  free-streara  tur¬ 
bulence  affects  an  aerodynamic  boundary  layer.  Therefore,  a  prediction  of  the 
turbulence  level  in  the  shear  layer  was  needed.  The  T'KE  model  was  changed 
from  one  based  on  an  average  value  of  TKE  across  the  mixing  zone  (Ref.  16)  to 
one  that  includes  the  solution  of  a  partial  differential  equation  for  the 
detailed  distribution  of  the  TKE  model.  The  modeling  follows  the  Prandtl 
Energy  Method  (Ref.  21). 

One  important  point  remained.  Since  the  parabolic  mixing  and 
burning  code  is  based  on  the  boundary  layer  form  of  the  equations  of  motion, 
no  upstream  influence  is  possible.  On  the  other  hand,  it  is  well  known  that 
heat  released  in  a  supersonic  flowfield  inside  a  duct  can  affect  the  initial 
conditions  (Ref.  10).  This  critical  behavior  was  included  by  using  the 
integral  analysis  of  Refs.  10  and  19  to  set  the  initial  conditions  and  p(z) 
along  the  duct. 

For  the  boundary  layer  calculations,  adoption  or  adaptation  of 
existing  codes  (or  parts  of  codes)  was  sought.  There  are  a  number  of 
elaborate  computer  codes  for  boundary  layer  flows  for  turbulent,  high-speed 
flow  conditions  (e.g.,  Refs.  22  and  23).  A  comprehensive  review  of  codes  for 
numerical  calculation  of  reacting,  laminar  boundary  layer  may  be  found  in  Ref. 
23,  but  none  is  suitable  for  a  scrarajet  combustor  because  the  chemistry 
generally  is  limited  to  air  (no  burned  fuel  species)  and  because  provisions  to 
account  for  the  rapidly  varying  boundary- layer  "edge"  conditions  under  turbu¬ 
lent  conditions  have  not  been  included.  Tl.e  necessary  extension  involved  not 
only  the  C-H-O-N  chemistry  but  also  significant  complications  in  terms  of  the 
viscosity,  Prandlt  number  and  Lewis  number  of  the  resulting  mixture.  There 
are  also  important  new  questions  with  regard  to  the  modeling  of  turbulent 
transport  processes  in  these  types  of  flows. 


The  modified  code  uses  the  basic  finite  difference  methods  of  Ref. 
22  for  the  numerical  solution  of  the  equations  of  motion  and  the  chemistry  and 
physical  properties  procedures  of  Ref.  25.  The  required  modifications  of  the 
code  of  Ref.  22  included:  a)  an  added  routine  for  the  solution  of  the  diffu¬ 
sion  equations,  b)  additions  to  the  routine  for  the  solution  of  the  energy 
equation  to  account  for  energy  transport  as  a  result  of  diffusion  and  for 
streamwise  variation  of  the  edge  conditions  on  total  enthalpy,  c)  revision 
and/or  replacement  of  routines  that  relate  enthalpy  to  temperature  and  species 
concentrations  and  those  that  calculate  the  equation  of  state,  d)  establish¬ 
ment  of  proper  calling  sequences  for  the  new  physical  property  subroutine,  and 
numerous  other  minor  modif ictions.  To  reduce  the  computing  cost  per  call,  the 
code  of  Ref.  25  was  substantially  modified  to  operate  as  a  subroutine,  and  the 
calculation  procedures  were  streamlined  and  simplified  to  the  C-H-O-N  system. 
The  routine  was  also  modified  so  that  the  transport  property  routines  could  be 
called  without  going  through  the  equilibrium  composition  calculations,  if 
desired. 


The  boundary  layer  flows  of  interest  differ  from  those  for  simpler 
aerodynamic  cases  in  two  important  ways.  First,  the  products  of  combustion 
are  diffusing  in  from  the  outer  edge.  This  changes  the  distribution  of 
density  and  physical  properties  across  the  whole  layer.  Second,  the  rapidly 
varying  boundary  layer  edge  conditions  result  in  profiles  of  atypical  shape  in 
the  outer  part  of  the  layer.  Hence,  it  is  inappropriate  to  apply  turbulent 
transport  models  that  have  been  developed  and  tested  only  for  conventional 
boundary  layer  flows.  A  model  more  fundamental  than  the  algebraic  eddy 
viscosity  or  mixing  model  of  Ref.  21  was  sought.  The  new  model  is  based  on 
the  variation  of  the  kinetic  energy  turbulent  fluctuations, 
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in  the  flow.  Prandtl  (Ref.  21)  suggested  the  form  of  one  such  approach.  A 
turbulent  viscosity,  y  ,  is  related  directly  to  k  with  a  length  scale,  l  ,  as 


The  variation  of  Z  is  specified  algebrica lly ,  but  a  separate  conservation 
equation  that  includes  k  must  now  be  solved  along  with  the  other  equations  of 
motion.  Importantly,  the  equation  for  k  involves  terras  for  the  production  and 
dissipation  of  turbulence  that  depend  upon  the  profile  shape.  Thus,  turbulent 
transport  is  made  sensitive  to  any  new  processes  that  influence  profile  shape 
(see,  e.g. ,  kefs.  26  and  27).  This  model  has  been  implemented  in  the  modified 
wall  boundary  layer  code. 

The  accurate  prediction  of  nozzle  flows  for  ramjets  and  scramjets  is 
important  for  efficient  design  and  accurate  system  performance  analysis.  The 
problem  is  made  complex  by  the  presence  of  strongly  nonuniform  nozzle  entrance 
profiles  typical  of  scramjets,  which  have  no  sonic  throat  at  the  end  of  the 
combustor.  Thus,  turbulent  mixing  and  perhaps  chemistry  can  be  expected  to  be 
underway  during  the  expansion  process.  The  rapid  expansion  common  in  nozzles 
of  practical  configuration  implies  substantial  radial  pressure  gradients 
also.  Taken  all  together,  these  requirements  imply  the  necessity  of  utilizing 
a  formulation  based  on  the  full  Mavier-S tokes  equations  for  a  compressible, 
multi-component,  turbulent  flow  with  finite  rate  chemistry. 

This  work  was  begun  with  the  VNAP2  code  developed  by  M.C,  Cline  at 
Los  Alamos  for  compressible,  single-component,  turbulent  flow  without  chem¬ 
istry.  Diffusion  equations  for  10  species  were  added,  as  were  additional 
terms  in  the  energy  equation  for  energy  transfer  due  to  mass  transfer.  Finite 
rate  chemistry  for  a  system  with  N-C-O-N  atoms  in  various  molecular  combina¬ 
tions  was  also  included. 

The  extension  of  the  VNAP2  code  to  consider  cases  with  varying 
composition  involves  a  considerable  effort.  For  example,  the  changes  to  the 
energy  equation  involve  not  only  the  additional  terms  for  energy  transfer  but 
also  additional  terms  for  varying  molecular  weight.  These  arise  because  the 
energy  equation  is  used  with  the  equation  of  state  to  develop  an  equation  that 
is  used  to  calculate  the  pressure. 

The  VNAP2  code  is  based  on  MacCormack's  explicit  method,  and  it  thus 
has  step  size  restrictions  from  the  flowfield  calculation.  Adding  provisions 
for  finite  rate  chemistry  introduces  step  size  restrictions  from  the 
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chemistry.  If  the  time  scales  involved  are  of  greatly  different  orders  of 

magnitude,  computational  instability  problems  can  result.  However,  from  a 
practical  point  of  view,  these  potential  problems  do  not  occur,  since  the  use 
of  this  code  can  be  restricted  to  cases  for  which  the  two  time  scales  are  of  a 
similar  order  of  magnitude.  If  the  two  time  scales  are  of  a  different  order 
of  magnitude,  we  have  the  following  alternatives.  If  the  chemical  times  are 
very  short  compared  to  flow  times,  we  assume  local  equilibrium.  If  the 
chemical  times  are  very  long  compared  to  flow  times,  we  use  finite  rate 

chemistry  calculations  by  themselves. 

The  Navier-S tokes  formulation  of  the  jet  raixing/combustion  code  is 
the  preferred  choice  at  this  time  because  it  includes  all  of  the  features  of 
the  parabolic  formulation  but  permits  pressure  feedback  and  radial  pressure 
gradients,  an  important  feature  when  computing  exit  nozzle  flowfields. 

Currently,  refinements  being  added  to  the  analyses  include  a  more  detailed 
description  of  the  base-like  flow  region  downstream  of  the  lip  of  the  gas 

generator  based  on  modifications  of  the  analyses  developed  in  Refs.  16-20  and 
28,  as  well  as  a  description  of  the  combustor/ inlet  interaction  region  once  it 
is  completed.  These  will  provide  more  detailed  initial  conditions  for  both 
the  boundary  layer  and  core  flow  models.  In  addition,  a  "two-fluid"  set  of 
equations  is  planned  for  incorporation  in  both  the  core  flow  and  boundary 
layer  codes  to  assess  the  effects  of  liquid  and/or  solid  particles  in  a  flow 
of  this  type.  Extension  of  the  f lamesheet/equilibrium  model  and  inclusion  of 
physical  property  routines  (to  avoid  the  binary  diffusion  approximation  cur¬ 
rently  used)  in  the  boundary  layer  code  are  also  planned. 

Resu 1 ts 

The  analysis  of  the  coaxial  mixing  and  combustion  region  is  used  to 
insure  an  adequate  length  of  combustor  for  complete  burning  and  to  predict  the 
level  of  non-uni f orral ty  of  the  combustor  exit  flow  entering  the  nozzle.  The 
primary  goals  of  the  wall  boundary  layer  calculations  are  predictions  of  the 
local  heat  transfer  and  skin  friction.  The  heat  transfer  must  be  known  to 
design  engine  -tooling  systems  and  structure,  and  the  skin  friction  in  the  com¬ 
bustor  wst  be  known  to  calculate  overall  system  performance  (e.g.,  gross 
The  r resent  wall  boundary  layer  analysis  Includes  the  important 


thrust). 


effects  of  entrainment  of  combustion  products  from  the  coaxial  flame  into  the 
boundary  layer. 
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Fig.  22  Predicted  combustor  wall  and  flame  sheet  contours  for  M0  =  4  and  7  flight 
with  ER  =  0.5. 


The  shape  and  length  of  the  two  flames  illustrated  in  Fig.  17  on  p. 
28  are  predicted  at  typical  Mach  4  and  7  flight  conditions  as  shown  in  Fig.  22 
on  p.  39.  The  longest  flame  is  for  CO  combustion;  it  is  about  1  meter  long, 
whereas  the  ^  flame  is  about  0.2  meters  long.  Fig.  23  on  p.  40  presents  the 
corresponding  predicted  central  core  temperature  and  species  profile  at  Mach  7 
with  a  fuel-air  equivalence  ratio,  ER,  of  0.5. 

The  edge  conditions  for  the  boundary  layer  calculations  were 
obtained  from  the  coaxial  jet  mixing  and  burning  code.  Fig.  24  on  p.  41  shows 
the  streamwise  variation  of  the  skin  friction  coefficient  and  wall  heat  flux 
for  Mach  7  flight  with  ER  *  0.5.  Here,  calculated  values  of  wall  skin 
friction  and  heat  transfer  show  that  both  increase  with  increasing  combustor 
heat  release,  a  result  observed  in  past  experiments  on  supersonic  combustion 
ramjet  engine  combustors  (Ref.  29).  The  sensitivity  of  the  local  heat 
transfer  to  the  local  character  of  the  pressure  distribution  is  also 
illustrated  by  the  steep  rise  of  the  local  heat  transfer  in  the  air  duct  ahead 
of  the  combustor  where  a  strong  shock  structure  is  located.  Fig.  25  on  p.  42 
presents  the  boundary  layer  CC>2  concentration,  temperature  and  velocity 
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profiles  at  the  combustor  exit. 
Refs.  15,  30  and  3L. 


Further  details  of  this  effort  are  given  in 
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Fig.  24  Combustor  wall  skin  friction  coefficient  and  heat  flux  distributions  for 
Mfl  =  7  flight  for  ER  =  0.5. 


Calculations  have  been  made  for  a  typical  scramjet  nozzle  case.  The 
nozzle  contour  chosen  for  analysis  is  plotted  in  Fig.  26.  The  inlet  flow  to 
the  nozzle  was  taken  as  the  outlet  flow  from  the  mixing  and  burning  code  for 
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Fig.  27  Nozzle  inlet  profiles. 


the  combustor.  i'he  inlet  property  and  species  profiles  are  shown  in  Fig.  27. 
An  "equivalent"  all-air,  uniform  inlet  flow  case  was  constructed  with  the 
outer  velocity  and  pressure  from  the  nonuniforra  case  and  the  same  total  mass 
flow.  Some  nozzle  property  exit  profiles  for  the  nonuniform  and  the 
associated  uniform  inlet  profile  cases  are  shown  in  Figs.  2S  and  29  on  p.  44- 
45,  respectively. 
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Fig.  29  Exit  profiles  -  uniform  inflow. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Combus tor/ Inle t  Interaction  Tests 

A  procedure  has  been  defined  and  tested  for  computing  the  instream 
flow  properties  in  a  ducted,  supersonic  flow  where  viscous  effects  are  impor¬ 
tant.  The  only  instream  data  required  as  input  in  this  algorithm  are  the 
pitot  pressures.  The  results  are  consistent  with  both  theory  and  experimental 
data  for  similar  flowfields.  This  consistency  is  with  respect  to  the  magni¬ 
tude  of  both  the  mean  (Ref.  2)  and  turbulent  fluctuation  components  (Ref.  3) 
of  the  various  flowfield  properties  of  interest.  Furthermore,  this  effort  has 
shown  that  severa-  alternate  algorithms  that  model  whe  instream  static  pres¬ 
sure,  rather  than  the  total  pressure,  were  not  as  successful.  Thus,  the  in- 
stream  total  pressure  p^(r,  z)  is  the  best  variable  to  model,  and  the  model 
presented  here  for  it  is  consistent  with  all  pertinent  physical  constraints. 

The  experimental  generation  of  the  simulated  precombustion  compres¬ 
sion  field  over  the  entire  range  of  test  variables  has  been  demonstrated.  The 
ducted  compression  field  is  generated  by  the  expansion  of  compressed  air  to 
pressures  that  are  considerably  less  than  the  ambient  exhaust  pressure.  By 
varying  the  annulus  plenum  pressure  over  a  narrow  range,  the  shock  structure 
in  the  annulus  can  be  moved  fore  or  aft,  with  the  accompanying  increase  or 
decrease  in  the  overall  static  pressure  rise  ratio.  The  generation  of  com¬ 
pression  fields  in  the  annular  isolator  duct  with  an  ambient  exhaust  and  with 
an  overall  pressure  rise  ratio  greater  than  the  normal  shock  value  has  been 
demonstrated  at  several  different  operating  conditions. 

Finally,  a  highly  efficient  interactive  computer  code  has  been 
developed  to  facilitate  experimental  data  reduction.  The  computer  graphics 
and  other  analytical  features  of  this  custom  designed  software  permit  one 
analyst  to  keep  pace  with  a  test  schedule  of  up  to  15  average  runs  per  week. 
An  average  run  involves  the  measurement  and  recording  of  about  200,000  data 
samples  (and  the  corresponding  200,000  time  samples)  for  about  250  locations 
in  the  flow.  Static  and/or  total  pressure  and  temperature  are  the  data 
measured. 
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The  experimental  investigation  is  now  armed  with  physical  insight, 
analytical  methods  and  data  reduction  codes  needed  for  rapid  and  significant 
progress.  The  difficult  period  of  grappling  with  unexpected  experimental 
data,  whose  implications  were  unclear,  is  over.  The  experimental  effort  is  in 
the  production  mode  of  operation,  where  the  direction  and  objectives  are  clear 
and  obtainable.  Both  the  hardware  and  the  software  are  fully  tested  and 
routinely  operational.  It  is  therefore  recommended  that  the  current  rapid 
progress  be  sustained  with  continued  funding. 

Component  and  Engine  Cycle  Modeling 

Computer  codes  have  been  developed  to  analyze  (1)  the  recirculation 
region  behind  the  gas  generator  base,  (2)  the  coaxial  mixing  and  burning  in 
the  supersonic  combustor,  (3)  the  combustor  wall  boundary  layer,  and  (4)  the 
nozzle  flow.  Base  pressure  predictions  for  Mach  2.03  and  2.30  frees tream 
conditions  have  been  completed.  These  calculations  have  included  the  effect 
of  the  upstream  growth  of  the  boundary  layer  momentum  thickness  on  the  base 
pressure.  The  coaxial  mixing  and  burning  code,  using  a  two  flame  sheet  model, 
has  predicted  that  the  H2  and  CO  flames  will  be  about  0.2  and  1.0  meters  long, 
respectively,  in  the  DCR  combustor.  Wail  shear  and  heat  transfer  distribu¬ 
tions  in  qualitative  agreement  with  observed  test  results  have  been  pre¬ 
dicted.  Finally,  the  numerical  modeling  of  the  ramjet  nozzle  flow  has  re¬ 

sulted  in  calculated  nozzle  exit  profiles  for  the  velocity,  density  and  Mach 
number  for  both  the  simplified  uniform  nozzle  inflow  condition  and  for  the 
real  nonuniform  entrance  condition.  The  nonuniform  entrance  condition  is 
specified  by  the  radial  variation  of  flowfield  properties  and  chemical  compo¬ 
sition  at  the  exit  of  the  supersonic  combustor,  as  predicted  by  the  mixing  and 
burning  code.  The  accurate  prediction  of  the  nozzle  flow  is  important  because 
the  net  thrust  developed  by  the  engine  is  highly  sensitive  to  the  magnitude  of 
the  nozzle  exit  stream  thrust  whose  value  in  turn  depends  on  the  finite  rate 
chemical  processes  in  the  nozzle. 

The  component  and  engine  cycle  modeling  effort  has  broken  much  new 
ground  but  has  not  yet  synthesized  all  the  modular  analytical  procedures.  The 
two  flame  sheet  coaxial  mixing  and  burning  code  needs  to  be  extended  to  treat 
multiple  species.  The  inclusion  of  both  finite  rate  and  equilibrium  chemistry 
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iato  the  mixing  and  burning  code  needs  to  be  completed.  The  preliminary  anal¬ 
yses  for  the  combustor  wall  shear  and  heat  transfer  need  to  be  improved  by 
replacing  the  initial  simplifying  assumptions  with  more  sophisticated  tech¬ 
niques.  The  analysis  of  the  recirculation  region  of  the  gas  generator  base,  a 
very  complex  and  challenging  problem,  requires  more  analysis.  The  effect  of 
the  lip  shock  that  follows  the  centered  corner  expansion  could  be  included. 

The  base  flow  analysis  is  very  important  because  it  produces  the  input  data 

for  the  mixing  and  burning  code,  which  produces  the  input  data  for  all  sub¬ 
sequent  analytical  modules.  The  extension  of  the  VNAP2  code  (from  Los  Alamos) 
for  the  numerical  simulation  of  ramjet  nozzle  flows  to  handle  varying  chemical 
composition,  to  include  additional  species  and  to  incorporate  new  procedures 
for  calculating  thermodynamic  quantities  should  be  completed.  It  is  therefore 
recommended  that  the  modeling  program  be  continued  in  order  to  permit  the 
synthesis  and  subsequent  productive  use  of  the  various  analytical  tools 
thathave  already  been  developed  and  to  permit  the  continued  refinement  and 
improvement  of  these  tools.  In  the  final  analysis,  these  recommendations  will 
prove  to  be  the  least  expensive  path  to  fully  develop  a  fundamental 

understanding  and  working  knowledge  of  the  flow  processes  in  supersonic 

combustion  ramjet  engines. 
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NOMENCLATURE 


A  =  area 

C£  *  skin  friction  coefficient 

ER  =  fuel-air  equivalence  ratio 

h  =  base  height 

k  =  kinetic  energy  turbulence 

l  *  length  scale 

M  *  Mach  number 

p  *  pressure 

=  wall  heat  flux 
w 

r  =  radius 

R  *  body  radius 

Re  =  Reynolds  number 

T  *  temperature 

u,  v,  w  *  velocity  components 

w  *  mass  flow  rate 

x,  2  =  axial  distance 

y  =  distance  from  wall 

Y  =*  mass 

p  *  density 

<5  ■  boundary  layer  thickness 

5*  *  boundary  layer  displacement  thickness 

6ke  *  boundary  layer  kinetic  energy  thicknes 

0  *  boundary  layer  momentum  thickness 

UT  »  turbulent  viscosity 

ip  *  stream  function 


*  supersonic  annulus 

■  edge  condition 

*  gas  generator 

*  jet 

■  maximum 

■  free  stream  (annulus) 

*  to  ta  1 

*  boundary  layer  momentum  thickness 


